Colossal non-saturating linear magnetoresistance in two-dimensional electron systems 

at a GaAs/AlGaAs heterointerface 
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Engineering devices with a large electrical response to magnetic field is of fundamental importance 
for a range of applications such as magnetic field sensing and magnetic read-heads. We show that a 
colossal non-saturating linear magnetoresistance (NLMR) arises in two-dimensional electron systems 
hosted in a GaAs/AlGaAs heterostructure in the strongly insulating regime. When operated at high 
source-drain bias, the magnetoresistance of our devices increases almost linearly with magnetic field 
reaching nearly 10,000% at 8 Tesla, thus surpassing many known non-magnetic materials that exhibit 
giant NLMR. The temperature dependence and mobility analysis indicate that the NLMR has a 
purely classical origin, driven by nanoscale inhomogeneities. A large NLMR combined with small 
device dimensions makes these systems a new and attractive candidate for on-chip magnetic field 
sensing. 



In a non-magnetic semiconductor, giant positive lin- 
ear non-saturating magnetoresistance (NLMR) can have 
both classical [H-Q and quantum origin [E], Q . The quan- 
tum NLMR, originally proposed by Abrikosov 0, Q , is 
applicable in the extreme quantum limit where Huj c 3> 
-Ef,/cbT [uj c and Ef are the cyclotron frequency and 
Fermi energy, respectively). This criterion can be at- 
tained in a restricted class of materials which include 
semimetals such as bismuth, or narrow band gap semi- 
conductors with very low effective mass (e.g. InSb |3(, 
graphene |7[ , and topological insultors 8] ) . The classical 
NLMR, on the other hand, is commonly observed in sys- 
tems with an inhomogeneous carrier (and hence mobility) 
distribution 0, It is a purely geometric effect, where, 
in the presence of transverse magnetic field, a misalign- 
ment between current paths and the externally applied 
bias mixes the off-diagonal components of the magne- 
toresistivity tensor, resulting in a NLMR. Several classi- 
cally inhomogeneous conductors, most notably the mildly 
doped silver chalcogenides (Ag 2 +aSe or Ag 2 +aTe) [l| and 
InSb polycrystals 3], display extremely large NLMR, 
where the inhomogeneity is associated with intrinsic dis- 
order such as grain boundaries, dopant clustering etc. 
Thus a handle on the disorder, both in magnitude and 
length scale, could yield a new class of high sensitivity 
magnetoresistive devices. However, achieving such a con- 
trol in bulk materials is not a trivial task. 

Two-dimensional electron systems (2DESs) in 
semiconductor multi-layers, in particular doped 
GaAs/AlGaAs heterostructures, offer a material 
platform in which disorder can be tuned with great 
precision using molecular beam epitaxy. At high carrier 
density, since the hetero-interface is physically separated 
from the ionized dopants, it is a homogeneous medium 



which hosts high mobility electrons. Therefore, it is not 
expected to be a good candidate for exhibiting NLMR. 
This is confirmed by numerous magnetoresistance (MR) 
measurements in 2DES, which are best studied in two 
limits. At high carriers densities, the MR of a 2DES is 
oscillatory in magnetic field (B) due to the Shubhnikov- 
de Haas effect, while at lower carrier densities, at the 
onset of localization, the MR increases exponentially 
with B due to variable range hopping (loj . Clearly, 
studies in neither regime have thus far revealed a NLMR. 

A 2DES however does become inhomogenous when it 
is depleted with a strong negative gate voltage (V g ) fllj | . 
This constitutes the backbone of our experiments in the 
following way: for a typical doped GaAs/AlGaAs het- 
erostructure with spacer thickness d, the Coulomb po- 
tential from the randomly scattered ionized dopants form 
the dominant component of disorder. This causes the 
conduction band minimum to fluctuate as a function 
of position (see Fig. lb) with a correlation length of 
~ 2d 12Hl4| . In this regime, the 2DES disintegrates into 



small puddles of charge that often manifest in Coulomb 
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blockade effects in mesoscopic devices [111 
essence, this increase in inhomogeneity arises due to the 
weakening of electrostatic screening of the background 
disorder potential landscape, causing the spatial fluctu- 
ations in carrier density to be of the same order as the 
carrier density itself. We show that such a strong density 
variation at the GaAs/AlGaAs interface does in fact give 
rise to a giant NLMR in a simple gate-tunable mesoscopic 
system. 

Our experiments are carried out in strongly insulating 
2DES's in the GaAs/AlGaAs heterostructure where we 
explore the magnetic response of the system as a func- 
tion of gate voltage (V g ) and drain-source bias (V s d). The 
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FIG. 1: (a) Schematic of the device — the gate modulates the electron number density only in region of 2DES below it 
(highlighted red) thus defining the mesoscopic region of interest (b) Top: Optical micrograph of the device. Scale bar is 20 /im. 
Bottom: Potential landscape profile for the electrons at very low number density (c) Inset: Equilibrium conductance as a 
function of the gate voltage. Region I represents gate voltage just above the pinch-off (the point when equilibrium conductance 
becomes zero); Region II represents another set of gate voltages below the pinch-off point. The main figure shows I — V 
characteristics of the device for these two distinct regions at 0.3 K. The arrow indicates the threshold voltage Vt- 



heterostructures contain a <5-doped layer of Si dopants of 
doping density 2.5 x 10 12 cm -2 placed 40 nm above the 
GaAs/AlGaAs interface. The 2DES lies 90 nm below the 



( a )io 2 



10' 



n 1 r~ 



2 I 

exp aB / 




FIG. 2: (a) MR of the sample for gate voltages above the 
pinch-off voltage (Region I) obtained using ac lock-in mea- 
surements (b) MR at different V s d below the pinch-off point 
(Region II). As source-drain is increased, the MR gradually 
transforms from a rapidly rising curve a to power law with 
exponent 1.1. The dashed line is R oc B. Inset: Percentage 
change in resistance AR/R (%) for two V s d of 14 mV and 
70 mV 



surface, with an as-grown mobility ~ 300 m 2 /Vs. The 
etched mesa and a surface gate dehne the effective ge- 
ometry of the device. Here we present results from a 
3 /imx7 [im device (see Fig. la and Fig. lb), though 
similar devices on two separate wafers showed qualita- 
tively similar results (see Supplementary Material). The 
gate is used to pinch the device off by application of a 
negative voltage, and the resulting conductance curve is 
shown in the inset of Fig. lc. Fig. lc shows / — V s d char- 
acteristics at two values of V g , representative of the two 
regions identified as Regions I and II. Region I represents 
the onset of localization at linear conductance G ~ e 2 /h, 
and has a reasonably ohmic / — V s d- Region II however 
lies deep in the pinched-off regime where the 2DES is 
expected to become inhomogeneous. Here, we find the 
/ — V s d characteristics are strongly non-linear with an in- 
sulating regime up to a threshold voltage Vt (indicated 
by an arrow in Fig. lc), which increases monotonically 
as V g is made increasingly negative (see Supplementary 
Material). Within this threshold, the current is immea- 
surably small (< 10~ n A), and conduction sets in rapidly 
only after the source-drain bias is increased beyond Vt- 
This rapid rise of current above threshold cannot be as- 
sociated with avalanche breakdown common in semicon- 
ductors IB- HI . This is because the applied electric fields 
are 3 orders of magnitude lower than the breakdown field 
in GaAs |T7( ■ The strong dependence of / — V s d curves on 



magnetic field (to be discussed later) also eliminates self- 
heating effects as a cause for the rapidly rising current. 
In fact, a systematic analysis of the / — V s d characteris- 
tics reveals that they follow a power law, characteristic of 
a disordered array of charge puddles (see Supplementary 
Material). This provides direct evidence that in Region 
II the system is indeed highly inhomogeneous. 

This inhomogeneity has a dramatic effect on the MR 
when subjected to a transverse magnetic field. At the 
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FIG. 3: (a) 7 — V s d characteristics as a function of magnetic field from T to 8 T (in steps of 1 T) at 0.3 K (Region II) (b) 
AR/R (%) for different V a d chosen from the shaded region of (a) (V s d step size is 5 mV) (c) A schematic illustrating the possible 
origin of NLMR: The light (dark) regions represent areas of high (low) carrier number density. White arrows represent the 
current paths. Left: At low Vsd, there is only one or very few current channels. Middle: For high V s d, the number of current 
paths proliferate rapidly. Right: On applying a perpendicular magnetic field at high V 3 d, the plurality of the paths in the 
disordered medium gives rise to non-trivial current trajectories which have significant transverse component (highlighted by 
dotted yellow arrows) that mixes the Hall contribution into the longitudinal MR, thus giving rise to NLMR (d) The mobility 
(calculated using Drude's formula) as a function of the conductivity of the device at 0.3 K, 5 K and 10 K. 



onset of localization (prior to pinch-off), represented by 
Region I in Fig. lc, the resistance R (obtained using equi- 
librium measurements) increases exponentially with B. 
As shown in Fig. 2a, R ~ exp(i? 2 ) at low B, which 
changes to R ~ exp(B) as B exceeds 2 T. These are 
characteristic features of hopping transport in the per- 
turbative (low-B) and non-perturbative (high-i?) regime, 
and recently studied in detail by some of us In 
the sub pinch-off regime, represented by Region II, R 
was found to behave very differently. Most notably, its 
structure is highly sensitive to the magnitude of V s d- 
To elucidate this, we define the dc magnetoresistancc 
R = R{V s d,B) — V s d/ 1 (V s d, B), and plot it as a func- 
tion of B for different values of V s d in Fig. 2b (we car- 
ried out the same analysis with differential resistance 
dV s d/dI(V s d, B), which did not yield significantly dif- 



ferent result). Strikingly, the variation in R weakens 
with increasing V s d, and for sufficiently high V s d » Vr 
[V T 5 mV at V g = -0.174 V, B = T), R increases 
linearly with B. In order to quantify the change in R, we 
evaluate the percentage change 

^(%). 1^4^^x100% 
R R( Vsd, 0) 

This quantity has been plotted for the two extremal V s d 
in the inset of Fig. 2b, which clearly shows that even the 
percentage change in R has transformed from a rapidly 
rising exponential curve to a strikingly linear form. 

To study this more systematically, we look at the evo- 
lution of the I — Vsd characteristics as B is increased from 
to 8 T. We have chosen a gate voltage corresponding 
to Vr ~ 30 mV (see Supplementary Material). Fig. 3a 
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shows that a non-zero B suppresses the current drasti- 
cally, leading to a positive MR that increases monoton- 
ically with B. Within the experimentally achievable B 
(8 T) the measured R did not show any sign of satura- 
tion. We have evaluated AR/R (%) (as defined above) 
for a few values of V s d above Vt (the range is highlighted 
in Fig. 3a). As shown in Fig. 3b, the MR in our meso- 
scopic 2DES reaches almost 10,000% at B = 8 T for 
V sd = 80 mV. As expected from Fig. 2b, the MR be- 
comes nearly linear with increasing V s d- The dashed line 
in Fig. 3b represents AR/R oc B 12 . 

The inherent inhomogeneity in charge distributions 
and a parabolic dispersion relation for the carriers makes 
the classical model for NLMR proposed by Parish and 
Littlewood 0, 0] particularly applicable in our devices. 
The findings in Ref. 4 followed from a numerical simula- 
tion of an equivalent node-link network model. Using the 
same physical principles, we have presented in the Sup- 
plementary Material an alternate theoretical analysis for 
NLMR in an inhomogeneous conductor, which augments 
existing theoretical descriptions [J, Isj ll9M22j . The clas- 
sical model requires current flow from source to drain 
to occur via multiple channels in order to realize the 
non-trivial magnetic response of NLMR. The process by 
which NLMR arises in our devices is depicted schemat- 
ically in Fig. 3c. At low V s d (left schematic), there are 
very few electron channels for conduction which are not 
sufficient in number to give rise to NLMR. However, a 
high V s d (middle schematic) opens up many more con- 
duction channels. This is similar to the non-equilibrium 
transport in disordered array of quantum dots where mul- 
tiplication of paths is directly connected to the conduc- 
tion threshold [23|. A perpendicular magnetic field dis- 
torts these current paths, which follows non-trivial tra- 
jectories through the inhomogeneous medium [4j , result- 
ing in a substantial transverse component (dotted ar- 
rows in right schematic). This allows a significant mix- 
ing of the the off-diagonal components in the magneto- 
resistivity tensor, thus leading to the NLMR. This quali- 
tative picture allows us to intuitively understand why we 
observe NLMR only at high values of V s d- We note that 
studies on mildly doped Si also reported a seemingly sim- 
ilar dependence of the quantity AR/R on source-drain 
bias [2j. However, the NLMR there was connected to 
an inhomogeneous electric field in the presence of space- 
charge injection. This scenario is certainly not applica- 
ble in our case since the bias applied in our experiments 
is significantly lower than that required to induce bulk 
semiconductor transport 171 Il8|. 

An important question that must be addressed is 
whether the inhomogeneity in our systems can quantita- 
tively explain the magnitude and field scales associated 
with the observed NLMR. Being a 5-doped heterostruc- 
ture, the dopant atoms are located at roughly the same 
distance from the 2DES, and hence the amplitude of con- 
duction band fluctuations is unlikely to vary widely from 
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FIG. 4: (a) AR/R at different temperatures: 0.3 K, 1.1 K, 
2.05 K, 5 K, 10 K for Region II. Dashed line represents 
AR/R oc B 12 (b) Comparison of MR of our device (high- 
lighted by dashed box) with other non-magnetic materials 
which have linear MR: Ag2+sSe [lj, Multilayer graphene Q, 
Bulk Si 0, InSb 0L LaSb 2 [H , Bi [H, HJ , Mg [23 , MnAs- 
GaAs composite [28f | . 



one location to another. Therefore, it is not unreason- 
able to assume that the width of mobility distribution 
Afi is smaller than the average mobility among the 
charge puddles. In such a case, the classical model pre- 
dicts a quadratic to linear MR transition at B ~ (a*) -1 , 
with AR{B)/R — c{fj,)B in the linear regime [22[ (c is 
a constant of order unity that depends on the details 
of mobility distribution). From Fig. 3(b), we find that 
AR(B)/R(0)B « 5 T -1 , which is consistent with the fact 
that the transition to linear MR occurs around 0.1 — 0.2 T 
(see Fig 2b). Moreover, a rough estimate of (/j) can be 
obtained from the mobility of the device at the onset of 
localization which is usually identified by conductivity 
— > e 2 /h for a 2DES. Below this conductivity, the system 
starts becoming disintegrated and incompressible [291 ] . In 
Fig. 3d we have shown the variation in the Drude mo- 
bility of our device as a function of its conductivity for 
a few temperatures. The limiting mobility at the on- 
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set of localization is about 5 m 2 /Vs, which is in good 
order-of-magnitude agreement with AR(B)/R(0)B. In- 
terestingly, (fi) is not a strong function of temperature in 
our heterostructures at least upto 10 K, presumably be- 
cause electron-phonon interactions contribute less to the 
resistivity than elastic residual scattering. This should 
also carry over to a weak temperature dependence of the 
NLMR. Fig. 4a shows that, indeed, the NLMR retains 
roughly the same magnitude for a temperature range 
from 0.3 K to 10 K, thus spanning nearly two orders 
of magnitude. 

Finally, we compare the sensitivity of the disordered 
GaAs/AlGaAs heterostructure with other non-magnetic 
systems that exhibit NLMR. The sensitivity of our device 
is extracted from its value at 1 T in Fig. 4a and is an 
increasing function of V s d (see Supplementary Material). 
The highest value of which that is directly observed is 
about 500% per Tesla at V sd = 70 mV. In Fig. 4b we have 
compiled the reported values of AR(B)/R(0)B (in per- 
centage per Tesla) from a wide variety of non-magnetic 
semiconductors and semimetals. Clearly, the sensitiv- 
ity of the GaAs/AlGaAs system is comparable to exist- 
ing non-magnetic materials. This, along with electrical 
tunability NLMR makes mesoscopic GaAs/AlGaAs het- 
erostructures unique in comparison to other materials. 
In fact, more efficient and clever gating assemblies could 
possibly achieve an even higher sensitivity and form a 
new class of magnetoresistive sensors. 
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Supplementary material: Colossal non-saturating linear magnetoresistance in 
two-dimensional electron systems at a GaAs/AlGaAs heterointerface 
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ANALYSIS OF THE I V sd CHARACTERISTICS 

Fig. la shows the gate voltage dependence of the / — V s d characteristics in Region II (as defined in the main text) 
at about 8 K. As is clear, the threshold Vr, which emerges in the / — V s d characteristics in the region II, has a distinct 
gate voltage dependence (Fig. lb). This is a remarkably linear dependence which does not break down up to the lowest 
gate voltages that we have applied. This linear dependence was also observed in other devices in a different wafer 
(not shown). The voltage threshold seems to be associated with the Coulomb blockade of the puddles of electrons, 
which is expected to form at such negative gate voltages. It is not clear how the linear dependence emerges, but may 
be the outcome the manner in which gate voltage changes the voltage barrier between these puddles of electrons. 
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FIG. 1: (a) Gate voltage dependence of the I — V s d characteristics of the device in Region II from V g = —170 mV to V g = —200 
mV in steps of 2 mV (b) Threshold voltage Vt as a function of gate voltage extracted from the / — V s d characteristics in (a) 

The fact that the system under study is highly inhomogeneous can be further ascertained by analysing the non- 
linear nature of / — V s d characteristics. In a typical / — V s d curve, the current rises rapidly after having crossed the 
threshold Vt- The I — V s d characteristics, in fact, display a power law with the form / = A(V s d — Vrr- Fig. 2a shows 
this for a few magnetic fields, at 4.9 K. Interestingly, very similar I — V s d curves have been observed in a disordered 
array of quantum dots [IH3- There, the bias threshold is directly related to the Coulomb gap of the series of quantum 
dots that define the electron trajectories from one lead to the other. £ depends on the type and strength of disorder, 
as well as the geometry of the device. It is not surprising that such / — V s d curves appear in our devices as well, since 
the 2DES is also expected to form a disordered array of electron puddles in the pinched-off regime at very negative 
gate voltages j^-Hoj. 

Interestingly, the value of £ drops sharply on the application of the magnetic field (Fig. 2b). At zero field, £ = 2.8, 
and falls to C — 1-8 by 1 T and saturates thereafter. A finer scan in another cool-down (inset of Fig. 2b) confirms 
that this behavior is indeed reproducible. The high ( = 2.8 that we observe at B = T is close to values reported 
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FIG. 2: (a)/ — V s d characteristics (with threshold Vt subtracted) cast in the log-log plot for magnetic fields T, 2 T, 4 T and 
8 T at 4.9 K. Dashed lines represent power law fitting of the respective curve (b) Plot of £ as a function of magnetic field B for 
measurements presented in (a). Inset: A finer scan of £ vs B (at 0.3 K) for another cool-down 

for array of quantum dots with structural disorder whereas £ ~ 1.6 at higher B is consistent with Middleton- 
Wingreen model of a capacitively coupled quantum dot array with on-site charge offset disorder jl[ . The reason for 
this transition of Q with B is unclear, but may be connected to redistribution of the charge inhomogeneity due to 
wavefunction squeezing in the puddles. Nonetheless, the / — V s d traces indicate that at sub-pinch-off gate voltages 
the 2DES is highly inhomogeneous, and may be visualised as an array of quantum dots with large structural and size 
disorder. 



BIAS DEPENDENCE OF NLMR 
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FIG. 3: AR/R as a function of V for magnetic fields from 1 T to 8 T at 4.9 K 



The fact that £ drops significantly from 2.8 (at B = 0) to 1.8 (at B = 0, and thereafter) implies that beyond 
the voltage threshold, / increases much faster when the magnetic field is zero. In other words, with increasing V s d, 
R(B — 0) drops faster than R(B > 0). As a result, AR/R (as defined in the main text) is an increasing function of 
V s d- Fig. 3 shows the variation of MR with V s d for all the magnetic fields from 1 T to 8 T. 
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FIG. 4: (a) Comparison of linear MR of the devl (reported in the paper) with that of dev2 (b) dev3 (c) dev4. In all plots, the 
dashed line is a power law curve with exponent 1.2 



LINEAR MAGNETORESISTANCE IN OTHER DEVICES 



Fig. 4 presents AR/R (calculated in the same way as in the main text) for other devices dev2, dev3 and dev4 in 
another GaAs/AlGaAs heterostructure wafer. Fig. 4a compares one of the MR curves for devl (from Fig. 3b of the 
main text) with those of dev2. The form of the MR for dev2 is same as devl and follows a power law with exponent 
1.2 (dashed line), thus exhibiting a near linear dependence. Fig. 4b and Fig. 4c show MR from two devices which 
show a slightly different exponent. 



THEORETICAL FOUNDATION OF NLMR IN AN INHOMOGENEOUS MEDIUM 
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FIG. 5: Schematic describing the model used to understand the NLMR: Left: A two-dimensional conducting medium is 
dispersed inhomogeneously with another material. Right: An effective medium equivalent to the one on the left. 



In order to understand the linear magnetoresistance in an inhomogeneous medium, we adapt the 3D magne- 
toresistance analysis of Stachowiak [ll| (itself motivated by Herring's effective medium treatment [13]) to 2D 
systems. Consider a two-dimensional conducting material in the x — y plane (Fig. 5) with longitudinal and Hall 
conductivities an and a ±2 respectively in presence of a magnetic field H z in the z— direction. Suppose we make 
another 2D conductor by inhomogeneously dispersing this material. The local conductivity of this composite material 

li 12 \ ^ -yy e dgjjQ^g the corresponding effective conductivity of the composite by a = { <T l 1 I 12 
-CTi2 an J y -0-21 a n 

We would like to obtain a relation between the local and effective conductivity. The measured resistivity of the 
sample is given by p — anj{^\\ + ^12) • 



is a = 



4 



Consider a small circular region of radius R made of the conducting material and embedded in an effective 
medium formed of the composite. Let there be a uniform electric field E 0x along the direction far from this 
disc. The resulting current j everywhere satisfies the continuity equation, V • j = 0. Using Ohm's law, j = erE etc., 
together with the definition E = — V^, where ip is the electrostatic potential, it is easily seen that the Laplace 
equation V 2 ijj — is satisfied. We need the solution for ip subject to the following boundary conditions: (a) 
linir-^oo V"0 = — (Eq x ,0), (b) ip(r) is continuous at r — R and (c) the normal component of the current at r = i?, 
j • f, is continuous. 

Solving the Laplace equation with these boundary conditions, one obtains the following solution for the potential 
tpi inside the disc: 



ipi(x,y) = Cnx + Dny, 

3(Tii(2<7ii + an)E 0x 



(25-n + on) 2 + (5-12 - cri 2 ) 5 



35-n(5-i2 - ai 2 )E 0x 

(2(7ii + Cll)- + ((712 - Cl2r 

Now we impose self-consistency, i.e., spatially averaging the internal field obtained above should give us the macro- 
scopic electric fields prevailing in the effective medium. Thus imposing (Ei x ) — Eq x and (Ei y ) = Eq v = we have 



35ii(25n + (7n) 



(25-n + CTn) 2 + (5-12 - CT12) 2 
35-n(5-i2 - (7i 2 ) 



, (25-n + an) 2 + (5-12 - (J12) 2 
For classical ohmic conductors, 



1, (2) 
0. (3) 



1711 ~ f + (H Z /H Q ) 2 ' 

= o-o (g, /H ) 
ai2 1 + (H z /H ) 2 ' 

where H^ 1 = er/m is related to the mean free scattering time r for the electrons in the material, and 00 = ne 2 r/m 
is the usual Drude conductivity. When H z 3> H 0} C12 ~ <r (H /H z ) ^> an ~ ao(H /H z ) 2 . This gives us a clue for 
the origin of linear magnetoresistance: in the inhomogeneous medium, the effective longitudinal conductivity takes 
on a bit of the local transverse conductivity. 

Consider now Eq. [3] The solution is ((712} = 5"i2. Next in Eq. [21 consider the strong magnetic field case first where 
we expect an 3> an on account of the admixture of 012 in the effective longitudinal conductivity. This simplifies Eq. 
[Uto 




which gives us 



*n »y !<(*12 -a 12 ) 2 )~0(l/H z ). (4) 

Correspondingly, the effective longitudinal resistivity p — an/i&ii + crf 2 ) ~ 0(H Z ). For weak magnetic fields, the 
effective longitudinal conductivity will be close to an and the magnetoresistance will go as H 2 . 



* Electronic address: aamir@physics.iisc.ernet.in 
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